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Spherically symmetric simulations of stellar core collapse and post-bounce evolution are 
used to test the sensitivity of the supernova dynamics to different variations of the input 
physics. We consider a state-of-the-art description of the neutrino-nucleon interactions, pos- 
sible lepton-number changing neutrino reactions in the neutron star, and the potential impact 
of hydrodynamic mixing behind the supernova shock. 



> 

cn 
a^ 

^ I Improvements of the neutrino-nucleon interaction rates 

' Recently, spherically symmetric Newtonian ||l], |2| , and general relativistic |^ hydrodynamical 

, simulations of stellar core-collapse and the post-bounce evolution including a Boltzmann solver 

I for the neutrino transport have become possible. No supernova explosions were obtained. 

Q^' Having reached a new level of accuracy of the numerical treatment, it is a natural next 

Q . step to reconsider the input physics that enters the models, remove imponderabilities of their 

pj I description and also test possible alternatives. Within the set of so-called "standard" neutrino 

^ ' opacities |0] which have been widely used in supernova models so far, for example, the rates 

of charged-current interactions of electron neutrinos and antineutrinos with free nucleons 
as well as neutral-current scatterings of neutrinos off free nucleons were calculated with the 
^ ' assumption that the nucleons can be considered as isolated, infinitely massive particles at rest. 

■ It is, however, well known that these interaction rates are significantly changed when energy 

transfer between the nucleons and leptons ("recoil"), nucleon-nucleon correlations due to 
Fermi statistics and nuclear forces, and other so far disregarded effects like weak magnetism 
corrections, are adequately taken into account (see Ref. ^]). Effectively, at densities p > 
10^3 g/cm^, the opacities become considerably smaller than in the "standard" approximation. 
Correspondingly, neutrino diffusion through the nascent neutron star and the emission from 
the neutrinosphere are enhanced and higher neutrino luminosities must be expected. 

Figure |^ shows results of Newtonian simulations of iron core-collapse and the post-bounce 
evolution of a 15 star using Boltzmann neutrino transport One model was cal- 
culated with an improved description of neutrino-nucleon interactions, which includes the 
detailed reaction kinematics and nucleon phase-space blocking, nucleon-nucleon correlations 
in the dense medium |8|, weak magnetism and the possible quenching of the axial 
coupling in nuclear matter For a reference calculation the conventional "standard" 
opacities Q were used. In both simulations we also included nucleon-nucleon bremsstrahlung 
A^A^ ^ uvNN (N € {n,p} denotes a free nucleon), which is the dominant production reac- 
tion of fi and r neutrinos and antineutrinos in the denser regions of the newly formed neutron 
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Figure 1: Comparison of the post-bounce evolution of a model computed with "standard" neu- 
trino opacities (black lines) and a model which employs the improved description of neutrino- 
nucleon interactions (red lines). Panel a shows the luminosities of Ve (solid lines), (dashed 
lines), and f^, v^-, i^^, Vt (individually; dash-dotted lines) as functions of time. In Panel b 
the radial positions of the shock (dashed lines) and the gain radius (solid lines) of the two 
models are compared. 



As expected from the arguments above, the model calculated with the new implementation 
of the opacities shows enhanced luminosities for Vf, and, even more pronounced, for Pg. This is 
caused by the fact that for i/g and recoil and nucleon correlations have the strongest effect 
at high densities, while for Vf, the weak magnetism reduces the opacities also at moderate 
densities p < lO^^g/cm^ [^, where most of the neutrino emission is produced during the 
time of consideration. In contrast, for Vf. both "corrections" counteract and almost cancel 
each other |^]. The higher luminosities (see Fig. |l]a) and somewhat larger mean energies of 
Vf. and Vf, increase the neutrino heating in the gain region below the hydrodynamic shock. 
Consequently, the shock propagates to larger radii when compared with the reference model 
(see Fig. %>). 

In view of the appreciable differences between both models we consider the inclusion of 
a state-of-the-art description of neutrino-nucleon interactions as a mandatory step towards 
more realistic supernova simulations, although the improved opacities alone do not lead to a 
successful supernova explosion in spherical symmetry. 



Lepton-number changing neutrino reactions 

Lepton-number changing interactions are expected in various extensions of the particle-physics 
standard model, notably in R-parity violating models of supersymmetry; see Ref. [jl^ for 
an overview and current experimental limits. Such processes might significantly affect the 
supernova dynamics [14|. As a proxy for this class of interactions we implemented the reaction 
fe + z/g + which opens a channel for "internal deleptonization" , i.e. quickly running 

down electron lepton number without diffusion to the stellar surface. This effect leads to a 
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Figure 2: Radial profiles of the electron fraction (Panel a, solid lines), total lepton fraction 
(Panel a, dashed lines), temperature (Panel b, dashed lines) and specific entropy (Panel b, 
solid lines) 50 ms after core bounce. Results of the models with lepton-number changing 
reactions are plotted as red (ct/usm = 10"'^) and blue lines ((t/itsm = 10"''). For comparison 
quantities of the reference model (ct/ctsm = 0) are shown as black lines. 



corresponding increase of the temperature of the stellar plasma by conversion of degeneracy 
energy to thermal energy (Fig. Q). The relevant time scale is fixed by the effective strength 
of the lepton-number changing reaction. We parametrize it by ct/ctsM) which is the lepton- 
number changing cross section normalized to the standard-model one for V(, + N ^ + N ■ 
We have calculated one model with ct/itsm = 10~^ where the deleptonization time scale 
is too slow to be effective during infall, and an extreme case where (j/ctsm = 10~^ and 
where coherent enhancement by the scattering off heavy nuclei is allowed, I'e + A ^ Ue + A. 
Despite of sizeable differences in the physical conditions in the core of the nascent neutron star 
(Fig. ^) we find a remarkable insensitivity of the overall dynamical evolution of the models 
to the dramatic modifications of the microphysics. In both models the hydrodynamic shock 
reaches a maximum radius of about 300 km at ~ 200 ms after bounce. These values are 
very close to the reference calculation with a/asM = 0. This finding is probably explained 
by the accelerated deleptonization and heating being confined to the innermost ~ 20 km 
of the neutron star (Fig. ^), because only there the lepton-number changing reactions are 
fast enough and the optical depth for these processes is larger than unity. During the first 
few hundred milliseconds after bounce most of the neutrino luminosity, on the other hand, 
originates from regions beween 50 km and 100 km where the temperature is essentially the 
same in all models (Fig. |2|b). Continuing the simulations over a period of nearly one second, 
we find slightly enhanced neutrino luminosities and mean energies and thus heating in the gain 
region at later times after bounce, but not to a degree that the conditions for explosions |15| 
would become more favourable. 
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Figure 3: Panel a shows the effect of our "mixing" algorithm on the temporal evolution of the 
radial profile of the specific entropy. From bottom to top the times correspond to profiles with 
increasing values of the postshock entropy. For the same model the trajectories of different 
mass shells (with a spacing of 0.02 Mq) are plotted as black lines in Panel b. The red line 
traces the radial position of the supernova shock. Note that the model includes only z^g and 
Pg, and can be directly compared with the simulation published in Ref. Q. 



Successful explosions by convective processes ? 

Multi-dimensional hydrodynamic simulations, coupled with simplified treatments of the neu- 
trino physics, have shown that large-scale convective overturn behind the supernova shock 
can aid the neutrino-driven mechanism and cause an explosion even if models fail in spherical 



symmetry (see Ref. [16| for an overview and references to original work). However, because 
of the uncertainties and approximations in the description of the neutrino sector, the actual 
role of convection in supernova models is not yet clear. Self-consistent and multi-dimensional 
simulations with a sufficiently accurate treatment of the neutrino transport have still to be 
carried out. 

Convective overturn has the helpful influence of suppressing energy losses by the reemission 
of neutrinos, because neutrino-heated matter expands outward instead of being advected 
through the gain region into the cooling layer. In addition, part of the matter falling through 
the shock can still accrete onto the neutron star, which helps keeping up a high value of 
the accretion luminosity. Furthermore, the expanding, neutrino-heated matter raises the 
pressure immediately behind the shock and thus drives the shock farther out. We try to 
mimick the latter effect of a multi-dimensional fluid flow between the shock and the gain 
radius by artificially (and instantaneously) mixing layers of negative entropy gradient behind 
the shock in our spherically symmetric simulations. This leads to a flattening of the entropy 
profile. We achieve this behaviour by setting the entropy in such a region to a constant value 
s{t,r) = so(t) (see Fig. ^), which is calculated in each timestep from the requirement that 
the internal (and thus also the total) energy in the affected region is conserved. Of course, this 
algorithm cannot accout for all the above mentioned effects associated with multi-dimensional 
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fluid motions, but at least one can test the sensitivity of the post-bounce supernova evolution 
to such a manipulation. 

As a consequence of the mixing and the thus increased postshock pressure, the hydrody- 
namic shock — which in the absence of mixing reached a maximum radius of about 350 km 
before it started to recede again — is indeed pushed out to a radius beyond 600 km (Fig. ^a), 
presumably leading to a weak supernova explosion. We interpret this result as an interesting 
hint that successful neutrino-driven supernova explosions might be in reach when state-of- 
the-art neutrino physics and a Boltzmann treatment of the neutrino transport are eventually 
combined in self-consistent multi-dimensional simulations. 
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